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Abstract 



In this paper, we study the spreading properties of the solutions of an integro-differential 
equation of the form ut = J*u — u + f{u). We focus on equations with slowly decaying dispersal 
kernels J(x) which correspond to models of population dynamics with long-distance dispersal 
events. We prove that for kernels J which decrease to slower than any exponentially decaying 
function, the level sets of the solution u propagate with an infinite asymptotic speed. Moreover, 
we obtain lower and upper bounds for the position of any level set of u. These bounds allow us to 
estimate how the solution accelerates, depending on the kernel J: the slower the kernel decays, 
the faster the level sets propagate. Our results are in sharp contrast with most results on this 
type of equation, where the dispersal kernels are generally assumed to decrease exponentially 
fast, leading to finite propagation speeds. 

Keywords: integro-differential equation; slowly decaying kernel; accelerating fronts; monostable; 
long distance dispersal. 
AMS: 47G20, 45G10, 35B40. 

1 Introduction and main assumptions 

In this paper we study the large-time behavior of the solutions of integro-differential equations with 
slowly decaying dispersal kernels. Namely, we consider the Cauchy problem: 



We assume that the nonlinearity / is monostable and that the initial condition uq is compactly 
supported. 

The equation (jl.ip arises in population dynamics [21] where the unknown quantity u typi- 
cally stands for a population density. One of the most interesting features of this model, compared 



ut = J*u-u + f{u), t>0, X G M 
u{0,x) = Uo{x), X gM. 




where J(x) is the dispersal kernel and 




to reaction-diffusion equations, is that it can take rare long-distance dispersal events into account. 
Therefore, equation (jl.ip and other closely related equations have been used to explain some rapid 
propagation phenomena that could hardly be explained with reaction-diffusion models, at least 
with compactly supported initial conditions. A classical example is Reid's paradox of rapid plant 
migration [IH [T2l [27] which is usually explained using integro-differential equations with slowly 
decaying kernels or with reaction-diffusion equations with slowly decaying - and therefore noncom- 
pact - initial conditions |25]. As we shall see in this paper, the use of slowly decaying dispersal 
kernels is the key assumption that leads to qualitative behavior of the solution of (jl.ip very different 
from what is expected with reaction-diffusion equations. 

Let us make our assumptions more precise. We assume that the initial condition : M — )• [0, 1] 
is continuous, compactly supported and not identically equal to 0. 

The reaction term / : [0, 1] — ?■ M is of class and satisfies: 

/(O) = /(I) = 0, f{s) > for ah s G (0, 1), and /'(O) > 0. (1.2) 

A particular class of such reaction term is that of Fisher-KPP type [n\ I19j . For this class, the 
growth rate f{s)/s is maximal at s = 0. Furthermore, we assume that there exist 6 > 0, sq £ (0, 1) 
and M > such that 

f{s) > /'(0)s - Ms^+^ for ah s G [0, sq]. (1.3) 

This last assumption is readily satisfied if / is of class C^'^. 

We assume that the kernel J : M — )• M is a nonnegative even function of mass one and with 
finite first moment: 

JgC°(M), J>0, J{x) = J{-x), I J{x)dx = land [ \x\J{x)dx < oo. (1.4) 

Furthermore, we assume that J(x) is decreasing for all x > 0, J is a function for large x and 

J'{x) = o{J{x)) as l^l — > +00. (1.5) 

This last condition implies that J decays more slowly than any exponentially decaying functions 
as |x| — )• oo, in the sense that 

Vr/>0, 3 x^eM, J{x) in [x^,oo), (1.6) 

or, equivalently, J(x)e''l^l — )• oo as |x| — t- oo for all e > 0. We shall refer to functions J satisfying 
the above assumptions ()1.4p . ()1.5p as exponentially unbounded kernels. 

The assumption p.Sp is in contrast with the large mathematical literature on integro-differential 
equations [21 [HI [151 [281 |30l [31] as well as integro-difference equations [HI [22], where the dispersal 
kernels J are generally assumed to be exponentially bounded as |x| — )• oo, i.e. : 

3ry > such that / J(x)e''l^l < oo. (1.7) 
Jr 

In this "exponentially bounded case", it follows from the results in |30j that, under our assumptions 
on no and /, the solution of (jl.ip admits a finite spreading speed c*. Thus, for any ci, C2 with 



< ci < c* < C2 < oo the solution u to (jl.ip tends to zero uniformly in the region |x| > C2t, 
whereas it is bounded away from zero uniformly in the region |x| < cit for t large enough. Thus, 
the spreading properties of the solution of (jl.ip when J is exponentially bounded are quite similar to 
that of the solution of the reaction-diffusion equation ut = Uxx + f{u) with n(0, •) = uq [3l HI [T71 [T9] . 
The existence of such a finite spreading speed is also true for other integro-differential equations 
with exponentially bounded dispersal kernels [2} [T ^ [25 1 [50]. 

Let us come back to problem (jl.ip with an exponentially unbounded kernel J. In this case, it 
is known that equation (jl.ip does not admit any traveling wave solution with constant speed and 
constant (or periodic) profile |32j. Moreover, numerical results and formal analytic computations 
carried out for linear integro-difference equations [2D] and linear integro-differential equations [23] 
indicate that exponentially unbounded dispersal kernels lead to accelerating propagation phenom- 
ena and infinite spreading speeds. In this article, we prove rigorously such results for the solution 
u of (jl.ip when the kernel J is exponentially unbounded i.e. J satisfies (jl.4p and (jl.5p . 

Our approach is inspired from |18j . where it was shown for a reaction-diffusion equation ut = 
Uxx + f{u) that exponentially unbounded initial conditions lead to solutions which accelerate and 
have infinite spreading speed. Here, we get comparable results starting from compactly supported 
initial data and with exponentially unbounded dispersal kernels. However, the interpretation of 
our results as well as their proofs are very different from those in [18j [25] . These differences are 
mostly due to the nonlocal nature of the operator tt i— )• J * u — u, and to its lake of regularization 
properties. 

2 Main results 

Before stating our main results, we recall that from the maximum principle |3Uj [52] and from the 
assumptions on uq, the solution u of (jl.ip satisfies 

< u{t, x) <1 for aU t > and x G M. 

For any A G (0, 1) and t > we denote by 

Ex{t) = {x G R, u{t,x) = A}, 

the level set of u of value A at time t. For any subset A C (0, J(0)), we set 

J~^{A} = {x G M, J(x) G A}, 

the inverse image of yl by J. 

Our first result says that the level sets Ex{t) of all level values A G (0, 1) (namely, the time- 
dependent sets of real numbers x such that u{t,x) = A) move infinitely fast as t — t- oo. 

Theorem 1 Let u be the solution of (jl.ip with a continuous and compactly supported initial con- 
dition no : M — ?• [0, 1] (uq ^ 0). Assume that J is an exponentially unbounded kernel satisfying (jl.4p 
and (jl.Sp . Then, 

Vc > 0, min u{t, x) as t oo (2.8) 

\x\<ct 



and for any given A G (0, 1), there is a real number tx > such that Ex{t) is non-empty for all 
t > tx , and 

mm{Ex{t)n[0,+oo)} - max{Ex{t) H {-oo,0]} 
lim = lim = +00. (^-y) 

t^+oo t t->+oo t 

Our next result gives a "lower bound" for the level sets Ex (t) in terms of the behavior of J at 



oo. 



Theorem 2 Under the same asumptions as in TheoUl for any A G (0,1) and e £ (0, /'(O)) there 
exists Tx^e > tx such that 

yt > Tx,e, Exit) C {(0,e"(^'(°)-^)*]} . (2.10) 
In our next result, we will either assume: 

Hypothesis 1 An exponentially unbounded kernel J satisfies Hypothesis [7] if and only if there 
exists cr > such that \ J' (x)/ J{x)\ is nonincreasing for all x > a and there exists Eq G (0, 1) such 
that 



[ Jizf^dz <oo. (2.11) 



or 



Hypothesis 2 An exponentially unbounded kernel J satisfies Hypothesis if and only if 



J'{x) 



J{x) 



O ( ) as \x\-foo. (2.12) 



Under these additional assumptions on the kernel J, we are able to establish an "upper bound" 
for the level sets Ex{t). 

Theorem 3 Let u be the solution of (jl.ip with a continuous and compactly supported initial con- 
dition uo : M — )• [0, 1] (no ^ 0). Assume that J satisfies either Hyp. [I] or Hyp. [H Then, there exists 
p > /'(O) such that for any A G (0, 1) there is Tx > tx such that 

yt > Tx, Exit) C { [e-P\ J(0)] } . (2.13) 

Theorem [2] together with Theorem [3] provide an estimation of the position of the level sets Exit) 
for large time t. In particular the inclusions ()2.10p and (j2.13p mean that, for any A G (0, 1) and any 
element xxit) G Exit), we have 

min ( J^i (e-(-^'(°)-")*) n [0, +oo)) < \xxit)\ < max (j-^ (e"''*) n [0, +oo)) , (2.14) 

for large t. 



3 Case studies 



Let us apply the results of Sec. [5] to several examples of exponentially unbounded kernels: 

• Functions J which are logarithmically sublinear as |x| — t- oo, that is 

J{x) = Ce-"l^l/i'^(l^l) for large (3.15) 

with a > 0, C > 0; 

• Functions J which are logarithmically power- like and sublinear as |x| — )• oo, that is 

J{x) = Ce-'^l^'l" for large \x\, (3.16) 

with a G (0,1), /3, C > 0; 

• Functions J which decay algebraically as |x| — )• oo, that is 

J{x) = Cl^r" for large (3.17) 

with a > 2, C > 0. . 

First, if J satisfies ([LlD and (l3J5]l then J satisfies Hyp. [D (but not Hyp. [2]). Theorem [2] and [3] 
then imply that for any level value A G (0, 1) and any e > 0, it exists p > /'(O) such that every 
element x\{t) in the level set Ex{t) satisfies: 

tln(t) < |xA(t)| < h\n{t) for large t. (3.18) 
a a 

Now, if J satisfies (II. 4p and (13.160 then J satisfies Hyp. [T] (but not Hyp. [2]) and it follows 
from Theorems [2] and [3] that the positions of the level sets Ex{t) are asymptotically algebraic and 
superlinear as t — )• +oo, in the sense that for e > 0, there is p > /'(O) such that 



■^'^"^ t'/" < \xxit)\ < I^^Y" t^/" for large t, 



(3.19) 



where x\{t) is any element of the level set E\(t) (see Fig. [T]). 

Next, if J satisfies (|1.4p and decays algebraically for large x as in (j3.17p . then J satisfies both 
Hyp. m and [2] and it follows from Theorems [2] and [3] that the position of the level sets E\{t) move 
exponentially fast as i — )• +oo in the sense that, for any A € (0, 1) and e > 0, there \s p > /'(O) 
such that ^ 

L!^l^t < ln(|a;A(t)|) < - t for large i, (3.20) 

a a 

for any x\{t) in the level set E\{t). The profile of the solution u{t, x) of (jl.ip with an algebraically 
decreasing kernel is illustrated in Fig. [5] (a) . 

We mention that Cabre and Roquejoffre [8] just established comparable estimates for the level 
sets of the solutions u of equations of the type ut = Au + f{u), where / is concave, uq is compactly 
supported or monotone one-sided compactly supported, and the operator A is the generator of a 




Figure 1: Plain line: position xo.2(0 of the level set -E'o.2(*) of the solution of with f{u) = 
u{l—u), uo{x) = max((l — (x/10)^), 0) and the exponentially unbounded kernel J{x) = (1/4)6"^^. 
Observe that xo.2(t) remains bounded by t H> J"^(e-^'(°)*) = (t-ln(4))2 and J-i(e-(/'(o)+V2)i) = 
(3t/2 - ln(4))2 (dashed lines) for large t. 

Feller semi-group. A typical example is the fractional Laplacian A = —(—A)" with < a < 1: if 
u is smooth enough and decays slowly to at infinty, 

(-A)"n(x) = c« / ^M^L^dy, for all x G M, 

where Cq, is choosen such that the symbol of (—A)" is In this case, the asymptotic exponential 
spreading of the level sets also follows from the algebraic decay of the kernel Jq,(x) = Ixl^*-^"*"^") 
associated with the operator A. We can notice that it is not a particular case of Theorems [2] and [3] 
since the kernel Jq, is singular at x = 0. 

Lastly, let us consider the example of a function J satisfying (jl.4p . (jl.5p and such that \J' / J\ is 
not monotone as |x| — >• oo, e.g. 

J(x) = C\x — sin(x)|~" for large |x| (with a > 2). 

Then J does not satisfies Hyp. 1, but still satisfies Hyp. 2. Thus, we can apply Theorems [2] and [3] 
which lead to the same estimates as (|3.20p . 

In all above examples the positions of the level sets increase super-linearly with time. This 
illustrates the accelerating behavior of the solution of (jl.ip for exponentially unbounded kernels. 
Coming back to Fig. [2] (a), we indeed observe that the distance between level sets of the same level 
tends to increase with time when time growths as i„ = an with a > and when J is exponentially 
unbounded, whereas it remains constant in the exponentially bounded case (Fig. [2] (b)). Moreover, 
in Fig. [2] (a), the solution tends to flatten as t — )• oo i.e. the lower the level A, the faster the level 
sets Ex(t) propagate. In particular, this implies that the solution does not converge to a traveling 
wave solution. This is coherent with the fact that (jl.ip does not admit traveling wave solutions 
when the kernel J is exponentially unbounded |32j . 
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(a) Accelerated propagation 




(b) Finite speed of propagation 



Figure 2: The solution u{t,x) of problem (jl.lj) at successive times t = 0, 3, . . . ,30, with f{u) = 
u{l — u) and uq{x) = max((l — (3;/10)^),0) : (a) with an exponentially unbounded kernel J{x) = 
(1 + (b) with an exponentially bounded kernel J{x) = (l/2)e~l^l. 



4 Proofs of the Theorems 
4.1 Proof of Theorem [1] 

We begin with proving that for any t > 0, liminf n(t,x) = 0. Let us define v(t,x) = v{t,x)e^* for 

>±oo 

all {t,x) G [0,oo) X M, where r = sup (/(s)/s) > /'(O) > and v satisfies the following problem: 

se(o,i] 

vt = J *v - V, t > 0, X gR, (4 21) 

v{0,x) = uo{x), X gW. 

Then v verifies Vt = J * v — v + rv on (0, oo) x M and v{0,x) = uq{x). From the maximum 
principle [301 [32] , we get < u{t,x) < v{t,x) for all {t,x) £ (0,oo) x R. 

Moreover, since uq is compactly supported and the operator n i— )• J*u—u is Lipschitz-continuous 
on L°°(R) n L^(R), the Cauchy-Lipschitz theorem implies that the solution t i— )• v(t,-) of prob- 
lem (|i:2ni belongs to C^([0, oo), L°°(M) n L^{R)). Integrating (jirm) over M and using ([LlD, we get 
II^(*)IIli(R) = II'"o||l1(]r)- This implies that u{t) belongs to L^(M) for all t > 0. Since u{t,x) > for 



any t > and x G M, we have: 

liminf n(t,x) = and liminf = for each t > 0. (4.22) 

Let us now prove that u{t, 0) — t- 1 as t — t- oo. Let / : [0, 1] — t- M satisfies ()1.2p . ()1.3p and such that 
/ ^ / in [0, 1], f{s) < f'{0)s for all s G [0, 1] and / is a nonincreasing function in a neighborhood 
of 1. We denote by u the solution of the Cauchy problem (|l.ip with the nonlinearity /. From the 
maximum principle u > it on [0, oo) x R. 

Then, let us set geis) = /(s) — £s in [0, 1], where e £ (0, 1) is small enough such that ^^(O) > 0. 

Set 

Ae = sup{s > I 5£ > in (0, s]} < 1. 

One can choose e > small enough so that < on (Ag, 1]. From (|1.4p we know that there exists 
> such that 



De = [ ^ J{y)dy = 1 - e. 
Let v_ be the solution of the following Cauchy problem: 



dtv = De{Je *v-v) + gdv), t > 0, X £ 
v{0,x) = uo{x), X G M, 



(4.23) 



where is a compactly supported kernel defined by: 

/ Jiy)dy 

We have 

dtv = ( J X \[-A,,A,]) * v - v + fin) < J *v-v + f{v). 

The maximum principle implies that < w(t,x) < u{t,x) < u{t,x) (< 1) for all t > and x G M. 
From Theorem 3.2 in [23], we know that v propagates with a finite speed c* > i.e. for all 

ce (o,c*), 



In particular, we have: 

Since A^ — )• 1 as e — )• 0, we get 



sup \v{t, x) - Ael — ;> as t — ;> +00. (4.24) 

\x\<ct 



lim v{t,0) = Xe< lim u{t,0) < 1. 



u{t, 0) — ^ 1 as t — > 00. 
It then follows that for any A G (0, 1) there exists a time ^ such that 

u{t, 0) > X for all t>tx. (4.25) 

Since the functions x 1— )■ u{t, x) are continuous for all t > 0, one concludes from ()4.22p and ()4.25p 
that Ex{t) is a non-empty set for all t >t\. 



Let us now prove (j2.9p . From [23j, we know that the propagation speed c* is the minimal speed 
of travehng wave solutions of problem (I4.23p . This speed verifies [9l [26] : 

4 = mm |i (^De ^ Je{z)e'^'dz - 1 + /'(O)) |. (4.26) 

Let j4 > 0, A G (0, 1), t > t\ and x\{t) E E\{t). Since J is exponentially unbounded in the sense 
of (|1.6|) . c* — oo as e — )• 0. Let us choose e > small enough such that > A and c* > A. Then, 
it follows from (j4.24p that 1x^(^)1 > At for t large enough. Since this is true for any A> {) and any 
x\{t) G Ex{t) we get (gS]) and dM]). □ 

4.2 Proof of Theorem [2] 

This section is devoted to the proof of a lower bound for min{£^A(i)n(0, +oo)} (resp. — v[i&yi{Ex[t)r\ 
(— oo, 0)}). The proof is divided into three parts. We begin with showing that the solution of (jl.ip 
at time t = 1 is larger than some multiple of J. Then, we construct an appropriate subsolution 
of (jl.ip which enables us to prove the lower bound for small values of A. Lastly, we show that the 
lower bound remains true for any value of A G (0, 1). 

More precisely, let us fix A G (0, 1) and e G (0, f'{0)). We claim that 

SA(t)C J-i{(0,e-(^'(°)-^)*]} (4.27) 

for t large enough. 

Step 1: u{l,-) is bounded from below by a multiple of J( ) Let us define 

v{t, x) = (no(x) + t{J * uo){x))e~\ 
Then it is easy to see that f is a subsolution of the following linear Cauchy problem: 

( = J *u-u, t > 0, X £R, 
1 u{0,x) = uo{x), X G M. 

Indeed, v{0,x) = uo{x) and for all {t,x) G [0, oo) x M: 

vt — J*v + v = —te~* {J * {J * uo)) (x) < 0. 

Since u is also a subsolution of the equation (11. ip verified by u, we get: 

u{t,x) > v{t,x) for all (t,x) G [0, oo) x R. 

Moreover, v{l,x) = (no(x) + {J *uo){x))e'^ for all x G M which implies that is exists C G (0, 1) 
such that v{l,x) > CJ{x) for all x G M. Finally, 



u{l, x) > v{l, x) > CJ(x) for all x G M. 



(4.29) 



Step 2: Proof of (j4.27p for small values of A We recall that there exist 5 > 0, sq £ (0, 1) and 

M > such that /(s) > /'(0)s - Ms^+^ for ah s G [0, so]- 
Define > by 

pi = /'(0)-6/2. 



Since J satisfies ()1.4p and ()1.5p . we can choose > such that for all |x| > 

\j'{x)\ < (e' 12) X J(x), 



where e' > satisfies 



Let us set: 



J( 

JO 



z)zdz < e/2. 



K= inf CJ = CJ(^i)>0, si = min (so, k), 

(-Ci.Ci) 



and 



B = max 



M V 1 + 



,-5 

'1 ' TT' 



M'(l + 5)(/'(0)-e/2) 



V 



> 0. 



Let g be the function defined in [0, oo) by 

g{s) = s- Bs^+^. 

We observe that 

g{s) < for all s > si and g{s) < si for all s > 0. 
Moreover, let < S2 < si be such that g'{s2) = and 



A2 = 5'('52)= max g(s) = — ^ ((1 + 5)5) 

sG[0,so] 1 + 



-1/5 



Let S,o{t) > be such that: 



C J{io{t))eP'' = S2 for all t > 0. 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



We can notice that for all t > 0, (,o{t) > and that (,o{t) is continuous and increasing in i > 0, 
since J is continuous and decreasing in [0,oo). 
Then, let us define u as follows: 



u{t, x) 



g{CJ{x)eP^') 



for all |x| > ^o(0 



A2 = g{s2) = g{CJ{Ut)V'') for all \x\ < Ut) 



for all t > 0. 



(4.36) 



Observe that < CJ{x)eP^^ < CJ(^o(i))e^^* = S2 < si when > Co(i), whence u{t,x) > for ah 
t > and X G M. Let us check that n is a sub-solution of (jl.ip . Since J{x) is nonincreasing with 
respect to |x| and u satisfies (|4.29p . we have 



u{0,x) = CJ{x) - B{CJ{x)Y+'' < CJ{x) < u{l,x) for all |x| > Co(0) 
u{0,x) = A2 = g{CJ{Co{0))) < CJ{x) < n(l,x) for all \x\ < Co(0). 



(4.37) 



Then, let us check that n is a subsolution of the equation satisfied by u in the region where 
u < A2. Let {t,x) be any point in [0,oo) x M such that u{t,x) < A2.AS aheady emphasized, one 
has < CJ{x)ef^^ < si, whence CJ{x) < si and \x\ > .^1 from (j4.32p . Furthermore, 

< u{t, x) < CJ{x)eP^^ < si < So < 1. (4.38) 

Thus, since / satisfies (jl.Sp . we get 

f{u{t, x)) > /'(O) (cj{x)eP^' - S(CJ(x))i+^e''^(i+'^)*) - M {C J {x))^+^ eP^^^+^'^K (4.39) 

Let us now show that J * u — u> — (e/2)n, for ah t £ [0, 00) and |x| > .^o(*)- Let t G [0, 00) and 

X > ^t) > 0. 

Remember that J is decreasing on [0, +00) and tha g is decreasing in [0, S2] (which imphes that 
u{t,y) is nonincreasing with respect to |y|). Then 

{J * u){t,x) — u{t,x) = / J{x — y){u{t,y) — u{t,x))dy 



> I J{x-y){u{t,y) -u{t,x))dy 

\y\>x 



> / J{x -y){u{t,y) -u{t,x))dy + J{x - y){u{t,y) - u{t,x))dy 

J —00 Jx 

f'QO 

> / {J{x-y) + Jix + y))iu{t,y)-u{t,x))dy. 



(4.40) 



Observe that for all y > x. 



ry 

^ ^ uitjU) — uit^x) = I dxu{t,s)ds. 



Furthermore, for all t > and s > (,o{t) 

d^u{t,s) = CJ'{s)eP''g' {CJ{s)eP'*) = CJ'{s)eP'' (l - (1 + 6) B {C J {s)Y eP^^' 
Since s > ^o{t) > : 

\dMt,s)\ < {e'/2)CJ{s)eP^'g' {CJ{s)eP^') 

< {e'/2)CJ{s)eP^' (l - B {C J {s)Y eP'^') = {e'/2)u{t, s) 

< {e'/2)u{t,x). 

Finally, for all y > x > ^o(i) 

u{t, y) - uit, x) > -{e'/2){y - x)u{t, x). (4.41) 



Then, from equation ()4.40p and (j4.4ip . we get 

POO 

{J *u){t,x) - u{t,x) > / {Jix-y) + J{x + y)){u{t,y)-u{t,x))dy 

J X 

> -{e'/2){j J{y - x){y - x)dy + j J{x + y){y - x)dy ] u{t, x) 

y 5oo poo ^ \ 

> -{e'/2)iyj J[z)zdz + j J{2x + z)zdzj u{t,x) 



> \yj J{z)zdz] u{t,x) 

> -{e/2)u{t,x). 



(4.42) 

The same property holds for x < —(,o{t) by symmetry of J and u with respect to x. It follows 
from (|i30D . (|i33D . (|i39D and (i02]l that, for aU t > and |x| > ^o{t) 

Uf{t, x) — {J * u{t, x) — u{t, x)) — f{u{t, x)) 

< piCJ{x)eP'^ - + S)B{CJ{x))^+^eP'^^+^^^ + e/2u{t, x) 

< {Pi - /'(O) + e/2) CJ{x)eP^' + (M - B[p^{l + 5)- /'(O) + e/2]) {CJ{x))'+'e''^^'+''>' 

< {M - B6pi){CJ{x)y+^eP^^'^+^^^ 

< 0. 

(4.43) 

Let us now check that u is a subsolution of the equation satisfied by u, in the region where 
u = X2- Let (t, x) be any point in [0, oo) x M such that u{t, x) = A2. The same arguments as above 
imply that for any point {t,x) satisfying u{t^x) = A2, i.e. for all {t,x) G [0, 00) x [— Co(0; Co(i)]; we 
get 



{J *u){t,x) -u{t,x) = / J{x-y){u{t,y) - \2)dy 

J{x- y){u{t,y) - \2)dy 



\y[>£.o{t) 

( J(x -y) + J{x + y)){u{t, y) - u{t, ^oit)))dy 



j POO POO \ 

> -{e'/2)[ J{y - x){y - Co{t))dy + J{x + y){y - Ut))dy ] X2 

/ POO POO \ 

> -{e'/2){j J{z + Coit) - x)zdz + J J {z + x + Coit)) zdz j X2 



> -e' J{z)zdz ) A2 

> -(e/2)A2. 



(4.44) 

For all X G [0,.^o(*)) it follows from (j4.36p that Ui{t^x) = 0. Let us show that this is also true when 
X = ^o(t). As already noticed, ^o(^) is an increasing function oft. Thus, for all /i > ^o(*) < ^0(^+^)5 



which imphes that u{t + h,^Q{t)) = A2 = M(i,^o(*))- As a consequence, 

h^O,h>0 h 



Moreover, 



j.^^ u{t,Ut))-li{t-h,Co{t)) _ X2-9{CJ{Co{t))eP'^'-'''^] 



h^O,h>0 h h^O,h>0 h 

For all > small enough and using the definition (|4.35|) of ^o(*)> we obtain: 

<7(CJ(eo(t))e^^(*~^)) = A2 - PiCJ{Ut))9\s2)h + o{h) = A2 + o{h). 
This implies that 

uAtMt)) = lim ^(^'^o(^))-f -^'^°(^)) = 0. (4.45) 

Since Ut{t,x) = for all x G [0,^o(*))) the above equality and the symmetry of the problem imply 
that: 

Ut{t, x) = for all |x| < ^o{t). (4.46) 
It follows from (lOD and that for all t > and |x| < Co{t), 

Ut{t, x) — {J * u{t, x) — u{t, x)) — f{u(t, x)) 

< (e/2)A2-(/'(0)A2-MA^+^) 

^ -w'(«)-/^)(i-(rTW(ra(rT^) j 

< 0. 

Using the above inequality together with (j4.37p and (j4.43p . the maximum principle implies 
that 

u{t -l,x) < u{t, x) for alH > 1 and x eR. (4.48) 
Fix now any real number uj small enough so that: 

<UJ < S2 

This real number uj does not depends on A but depends on e, as well as on J and /. Remember that 
> is such that E^^ (t) is a non-empty set for all t >tuj. Since J is continuous and decreasing on 
[0, +00), there exists then a time tu) > max {t^}, 1) such that for all t > t^, it exists ya;(i) £ (Ci, co) 
such that 

CJ(y.(t))e^^(*-i)=u.. 

Furthermore, the function yu]{t) : [to;, 00) — )• [^1,00) is increasing and continuous. 
Lastly, let be the open set defined by 

n = {{t, x) e (L, +00) X M, |x| < yUt)}- 



We claim that inf u > 0. Indeed, if {t,x) £ 0, is such that CJ{x)e''^^^~^^ > S2, then \x\ < ^o(i — 1) 
and 

u{t - = A2 = g{s2) > g{uj) > 0. 
Otherwise, (t,x) is such that oj < CJ{x)eP^^^^^^ < S2, whence |x| > ^o{t — 1) and 

u{t -l,x)= 5((7J(x)e^i(*i)) > g{uj) > 0. 

Finally equation (I4.48P implies that 

u{t,x) > g{uj) > for aU {t,x) G Q. (4.49) 

Thus, setting 6 = g{uj), we get that if A G (0,6*) and if x G Ex{t) for t > max(t;\j, t^), then 

|x| > yUt) > 6 > Co(0) 

Since pi = /'(O) — e/2 > /'(O) — e, there exists then a time Tx,e > max(tA,iaj) such that, 

Vt > Ta,,, Vx g ii;A(t), J(x) < J{yUt)) = ^e-''^* < e-^^'^"^)*. (4.50) 
This proves KTHi for A G (0,6*). 

Step 3: Proof of ()4.27p for any A G (0,1) Assume that A G (0,1). Let Ug g be the function 
defined by 

/ 0(1 -|x|) if |x| <1, 

Un nix) = < 

1 if |x| > 1, 

where 6 = g{uj) is given in Step 2. Let us set f{s) = f{s) — (1 — Dx)s for all s G [0, 1] with 

Dx= / J {y)dy e {0,1), 

where ^1 is chosen large enough such that /'(O) > and /(s) > in (0, A]. 
We consider the solution Ug of the Cauchy problem: 



j dtug = Dx{Jx*Ue-yie) + fiHe), t > 0, x G M, 
\ Me(0,x) =2^^,0(2;)' 2; G M, 

where Jx is a compactly supported kernel defined by: 

Mx) = l[_^,,^,](a;), X G M. 

/ J{y)dy 

It follows from (j09]) that 

VT>L, V|e| <yUT)-l, VxGM, u{T,x)>Ugn{x-0), 



(4.51) 



whence 

VT > L, V|^| < yu,{T) - 1, yt> 0, Vx G M, ^(T + x) > Ug{t, x - ^) (4.52) 
from the maximum principle. Indeed, we have for all (t, x) G (0, oo) x M: 

X) - {J * Ujg{t, X) - Ug{t, X)) - fiUffit, X)) 

< dtugit, x) - DxJx * ue{t^ x) + us{t, x) - f{ue{t, x)) 

< dtug{t,x) - Dx {Jx *U0{t,x) -Ujsit^x)) - f{ujg{t,x)) 

< 0. 

Moreover, we know from Theorem 3.2 in [23j that there exists > such that 
liminf inf x) = sup{s > | / > in (0,s)} > A. 

t->-oo \x\<clt 

In particular, there exists T_x > such that Ujq{Tx,x) > A. Therefore, ()4.52p implies that 

VT > L, V|x| < yUT) - 1, u{T + Tx, x) > A. 

As a consequence, there exists Ta e — ^^^(tu! + Tx,tx) such that for all t > T_xi, and for all 
x G Exit), one has \x\ > y^{t — Tx) — 1 and 

Jix) < J it - Tx) - 1) = ^^e"-(-^^) X 'f^l'^i''^ < e-^rio)-e)t^ 

since — Tx) — )• oo as t — )• oo and '^^ — 1 as s — )• oo, from (II. 5p . This implies (I4.27P and 

J(s) 

completes the proof of Theorem [2j □ 
4.3 Proof of Theorem H 

In this section, we prove an upper bound for max{£'A(i) 1^(05 +°*^)} (resp. — mm{Exit)rii—oo,0)}). 
The proof of this upper bound is based on the construction of suitable supersolutions of (jl.lh . The 
construction of such supersolutions strongly relies on Hypotheses [1] and [2j 
We shall prove that there exists p > such that, for any A G (0, 1), 

Exit) C J'^ {[e-P\ J(0)]} for large t. (4.54) 

Since / G C-'^([0, 1]), the "per capita growth rate" /(s)/s is bounded from above by r = 
sup ifis)/s) > 0. 

sg(0,l] 



Proof of (|4.54p under Hypothesis [T] Assume that J satisfies Hyp. [TJ Then J' / J is negative 
and nondecreasing on [a, +oo). Therefore ln(J) is a nonincreasing convex function on [a, +oo). It 
is then possible to define a function (/? : [0, oo) — [0, oo) and r G [0, a] such that (/? is nondecreasing 
and concave and 



(fiO) = 0, v3(+oo) = +00, and (fix) = (eo — 1) ln( /(x + r)) on [a — r, +oo). (4.55) 



By concavity, we have the following property, for all y > x > 0: 



ip{y) - ip{x) < ^p{y - x). (4.56) 

Thus, we claim that: 

Vx > 0, Vy G M, ip{x) - ip{\x - y\) < ip{\y\). (4.57) 

Even if it means increasing a, one can assume without loss of generality, that J < 1 on [ex, oo). 
Indeed, if y < then (^(x) — ip{\x — y|) < < (/^dyj). If y > 0, since f is nondecreasing and 
from (|4.56|) . we have ip{\y — x\) > (/?(max(x, y)) — (/?(min(x, y)) > (p{x) — (p{y). Notice that (|4.57|) 
implies immediately that 

Vx G R, Vy G M, (^(|x|) - ip{\x - y|) < (^(|y|). (4.58) 

Let us define 

Mx) = e-'^(l^l) for all x G M. 



Using (|4.58|) . we get: 



(J*0)(X) _ f H^-y\j^y^^y 



(t>{x) Jr <P{x) 

j(y)e'^(l^l)-v(l^-2^l)dy (4.59) 

< / Jiy)e^^\y\Uy= [ J{y)/mdy. 



Moreover, 



J{y)/<t^{y)dy= j J{y)e'P^\y\)dy = [ J{y)e'^''^y\Uy + [ J{y)e^^\y\Uy 

JR J\y\<(7 J|»/|>o- 

= / J{y)e'P^\y\^dy+ f (i^lllT 'j^yfo < 

J\y\<a J\y\>a V J iy) ) 

(4.60) 

J(y _|_ ■y\ 

from Hyp. [1] and since -— )• 1 as y — )• ±00 from (II. 5p . 

Since uq is compactly supported, there exists (Ti > such that ^{x) > uq{x), for all |x| > ai. 
Finally, set 

po = max I J J{y)e'^^\y\Uy, l| - 1 + r, 

and define u as follows: 

V(t, x) G [0, +00) X R, u(t, x) = min ( ^^e''«*, 1 ) . 

Observe that ^0(2;) ^ m(0, x) for all x G M. Let us now check that u is a supersolution of the 
equation (jl.ip satisfied by u. Since u < l,it is enough to check that ?I is a supersolution of (jl.ip 



whenever u < 1. Note that since (j){x) is nonincreasing with respect to u{t,x) < 1 imphes that 
\x\ > ai. Assume that {t,x) G [0, +00) x [o"i,oo) and u{t,x) < 1, then it foUows from (I4.59P that 

{J*u){t,x) < {J*(j)){x)- 



^ 6(x) ^ JL 



[ J{y)/mdy 



= u{t,x) / J{y)/(p{y)dy. 

This imphes that for all {t, x) such that u{t, x) < 1 

ut{t, x) — {J * u){t, x) + u{t, x) — f{u{t, x)) 

> pouit, x) — {J * u){t, x) + u{t, x) — ru{t, x) 

(4.61) 

> (Po- / J{y)/cPiy)dy + l-r)u{t,x) 

> 0, 

from the definition of po- The parabolic maximum principle 130\ I32j then implies that: 

u{t,x) < u{t,x) < ((/)(x)/(/>(f7i))e^"*, for ah {t,x) E [0,oo) x M. 

For all t >tx (so that Exit) is not empty) and for all x £ Ex{t), there holds 

(|x| < ai) or (|x| > fji and A = u{t,x) < {(l){x)/(j){ai))eP°^) . 

In all cases, we get that 

Vt > tx,yx G Ex{t),(t>{x) > min((/>(CTi), A(/)(o-i)e"''°*). 

Then, from the definition of (j) for large x, and since J(s + t)/J{s) — )• 1 as s — >• ±00 for any 
p> Pq/{1 — Eq), there exists a time Tx > tx such that 

Vt > Tx, Vx G Exit), J(x) > e-P\ (4.62) 

which gives ()4.54p . 

Proof of (|4.54p under Hypothesis [2] Assume that J is satisfies Hyp.O Since J(x) is decreas- 
ing with respect to |x|, we get the following inequality for any x > : 

iJ*J){x) Jix-y) r Jix-y)Jijj) Jiy) 

-JW- = L-Jixr'^'^''^Jo Ax) ''^L Jix)'^'-'^'' 

< 1+ r^'itzylj(y)dy+ r Ii^Jiy)dy (4.63) 

Jo 'Ji.x) Jx/2 J{X) 

Jix/2) 
- + J(x) • 



From the symmetry of J, the inequahty also holds for x < 0. Moreover, since J satisfies ()2.12p . 
there exists > 0, and Co > such that 



In 



j{x/2y 



= I —rrr ds < — ds = Co In (2) for all Ixl > xo- 

/ J\x\/2 J\x\/2 s 



This implies that there exists K > such that for all x G M, 

^-^fM<l + K. (4.64) 
J(x) 

Since uq is compactly supported, there exists cti > such that J(cri) < 1 and J(x) > uo{x), for all 
|x| > (Ti. Then, set po = r + K and for all {t, x) G [0, +oo) x R : 



u{t, x) = min 



Observe that ^0(2;) < n(0, x) for all x G M. Let us now check that n is a supersolution of the 
equation (ll.lj) satisfied by u. In the region {t,x) such that u{t,x) = 1, the same arguments as in 
Sec. 14.31 lead to: 



ut{t, x) — J * u{t, x) + u{t, x) — f{u{t, x)) > 1 — J * u{t, x) > 0. 

Let us check that u is also a supersolution of (jl.ip when n < 1. If t > 0, |x| > cji and u{t,x) < 1 
then it follows from (|4.64p that 

{J*u){t,x) < (J*J)(x) 



gPO* 



< {l + K)u{t,x). 



This implies that 

ut{t, x) — ( J * u){t, x) + x) — f{u(t, x)) > pou{t, x) — {J * u)(t, x) + u{t, x) — ru{t, x) 

> {po-{l + K) + l-r)u{t,x) 

> 0. 



(4.65) 

The parabolic maximum principle |30l [32] implies that: 

J(x) 

u{t,x) < u{t,x) < -T7—Te^°\ for all (t,x) G [0,oo) x M. 

For all t >tx (so that Ex{t) is not empty) and all x G Ex{t), there holds 

/ J(x) 
{\x\ < fJi) or |x| > (Ti and A = u{t,x) < e^"* 

V d{(^i) 

In all cases, one gets that 

yt > tA,Vx G Ex{t),J{x) > min(J(CTi),AJ(o-i)e-^«*), 
Then, for any p > po> /'(O) > 0, there exists a time Tx > such that 

> Ta, Vx G Ex{t), J(x) > e-P\ (4.66) 
which proves ()4.54p . □ 



5 Discussion 



We have analyzed the spreading properties of an integro-differential equation with exponentiahy 
unbounded or "fat-tailed" kernels. Since the pioneering work of Kot et al. [20], there have been 
few mathematical papers on integral equations with exponentially unbounded kernels. However, 
such slowly decaying kernels are highly relevant in the context of population dynamics with long 
distance dispersal events [HI [T2] . 

We proved that for kernels J which decrease to slower than any exponentially decaying function 
the level sets of the solution u of the problem (|l.ip propagate with an infinite asymptotic speed. 
This first result shows the qualitative difference between dispersal operators with exponentially 
unbounded kernels and dispersal operators with exponentially bounded kernel which are known to 
lead to finite spreading speed [21 [151 ISSl [30] . This result supports the use of "fat-tailed" dispersal 
kernels to model accelerating propagation or fast propagation phenomena [HI [121 [2Z] • 

Moreover, we obtained lower and upper bounds for the position of any level set of u. These 
bounds allowed us to estimate how the solution accelerates, depending on the kernel J: the slower 
the kernel decays, the faster the level sets propagate. Through several examples, we have seen in 
Sec. [2] that the level sets of the solution of problem (jl.ip move almost linearly when J is close 
to an exponentially bounded kernel (see example (j3.15p . J{x) = C'e"'^'/''^^!^'!-* for \x\ ^ 1) while 
the level sets move exponentially fast when the kernel J has a very fat-tail (see example (|3.17p . 
J{x) = C|x|-" for |x| > 1). 

It is noteworthy that our results have been derived under assumptions more general than the 
KPP assumption f[s) < f'{0)s. Indeed, results of Theorems [H [2] and [3] hold with nonlinearities / 
which may take a weak Allee effect into account: the maximum of the "per capita growth rate" 
f{s)/s is not necessarily reached at s = 0. In ecological models, the Allee effect can occur for 
various reasons [1]. For instance, at low densities individuals may have trouble finding mates. Our 
spreading properties in the case with a weak Allee effect are in agreement with the numerical results 
in [29j which show that exponentially unbounded dispersal kernels can lead to infinite spreading 
speeds. The conclusion is very different when the nonlinearity / takes a strong Allee effect into 
accounts, that is if /(s) < for small value of s. Indeed it is proved in this case that problem (jl.ip 
admits traveling wave solutions with constant speed [5l[10l[13]. Thus the solutions of problem (jl.ip 
with such nonlinearities have a finite speed of propagation. 

One could wonder whether the solution of problem (jl.ip with exponentially unbounded kernel 
converges to some kind of "accelerated traveling wave solution", that is a solution u{t,x) = (/)(x — 
c(t)) where c is a superlinear function. Our numerical computations suggest that the answer is no 
(see Fig. [2| (a)) since the solutions u becomes flat at large time. The computations also suggest 
that the solution u{t,x) is not a generalized transition wave in the sense of [H [7]. Such a result 
has been proved by Hamel and Roques [18] for the solution of a reaction-diffusion problem with an 
exponentially unbounded initial condition. In our case, the proof seems to be more involved since 
the operator u^J*u — u is not a differential operator. 
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